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Abstract. Three methods for studying stationary measures of stochastic differential 
equations with jumps are considered. These stationary measures are given by Markov 
measures, solutions of Fokker-Planck equations, and long time limits for the distributions 
. of system states. 
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1. Introduction 

Stationary measures for stochastic differential equations (SDEs) and invariant measures 
for Markov processes have been studied extensively, as in [2J El El HS1 ELY], among others. 
Moreover, there are interesting relationships between stationary measures and invariant 
measures, such as the one-to-one correspondence between the set of invariant Markov 
measures (c.f. Definition 13.11 in Section 3) and the set of stationary measures (c.f. [31 
|5]), as well as the correspondence between the solutions of Fokker-Planck equations and 
stationary measures (c.f.[13j). These results play an important role in the development 
CN ■ of the theory for random dynamical systems associated with SDEs. 

SDEs with jumps and random dynamical systems associated with them are considered 
by a number of authors (c.f. [I] [2] [Z] [ID] [TT] [J3] [IT]). We note that Albeverio- 
^ ' Riidiger-Wu pQ discussed stationary measures for SDEs with jumps, in the context of 

Levy type operators and considered mainly infinitesimal invariant measures. The concept 
of infinitesimal invariant measures is weaker than that of the usual stationary measures. 
Moreover, Zabczyk [17] studied stationary measures for linear SDEs with jumps. 

We ask, naturally, whether the above correspondences for usual SDEs also hold for 
SDEs with jumps. This question will be answered in Sections 3 and 4 in the present 
paper. 

Now we briefly sketch our method. We begin from SDEs with jumps as random dy- 
namical systems and examine their Markov measures. Then we consider SDEs which 
are driven by Brownian motions and by a-stable processes. By a functional analysis 
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technique, stationary measures are investigated. When the coefficients of the SDEs are 
sufficiently regular, the long time limits for the distributions of the solutions are shown 
to be stationary measures, with the help of Malliavin analysis. 

This paper is arranged as follows. In Section 2, we introduce random dynamical systems 
and related concepts. Sobolev spaces and a-stable processes are also introduced. The 
content to obtain stationary measures from Markov measures is in Section 3. In Section 
4, we consider special stochastic differential equations driven by Brownian motions and te- 
stable processes. In Section 5, we study SDEs with jumps, under certain regular conditions 
on the coefficients. 

The following convention will be used throughout the paper: C with or without indices 
will denote different positive constants (depending on the indices) whose values may 
change from one place to another. 



2. Preliminary 

2.1. Random dynamical systems and related concepts. Let D(M + ,M d ) be the set 
of all functions which are cadlag, defined on M + and taking values in M d . We take sample 
space Q = D(R + ,~R d ). This sample space becomes a complete and separable metric space 
(c.f. [5]), when endowed with the following Skorohod metric d: 



d(x,y) :-- 




A(t)-M.) 



m=l ) 



t-S 

for all x, y G O, where x m (t) := g m (t)x(t), y m (t) := g m (t)y(t) with 

1, if t ^ m, 
9m(t) ■= { m + 1 — t, ii m < t < m + I, 
0, if t ^ m + 1, 



and 



d° m (x,y) ■■= sup \x(t) - y(X(t))\ . 



Moreover A denotes the set of strictly increasing and continuous functions from M + to 
1R + , and a A b := min{a, b}. We identify a cadlag function u(t) with a (canonical) sample 
to in the sample space Q. 

The Borel cr-algebra in the sample space Q under the topology induced by the Skorohod 
metric d is denoted as T. Note that T = cr(u(t),t e IR + ) (c.f.[6j). Let P be the unique 
probability measure which makes a canonical process a Levy process, and denotes (J-"t) tG R + 
the complete natural filtration with respect to P. Thus (f2, J 7 , P; (•7 r t)te]R + ) is a complete 
filtered probability space. 

Define for each t £ R + , 

(9tw)(-) =u(t+ ■) -u(t), uett. 
Then 9 is a one-parameter semigroup on Q, Q is invariant with respect to 9, i.e. 

e^n = n, for ail t e r+, 

and P is ^-invariant, i.e. 

P(9; 1 (B)) = P(B), for all B e 7 ,t G R+. 
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Thus (ft, J 7 , P, (Ot)tm+) is a metric dynamical system (DS), as defined in [3J. This metric 
DS is in fact ergodic, i.e. all measurable invariant sets have probability either or 1. 
(c.f.0) 

Definition 2.1. A measurable random dynamical system on the measurable space (X, B) 
over a metric DS (Q, J 7 , P, (9t)teR+) with time R + is a mapping 

ip : R + x x X 4 X, (t, cu, x) \-> <p(t, u, x), 

<f(t, u) := (p(t, u,-) : Xh-X, 

such that 

(i) Measurability: is £>(M + ) £g> T <S> B / B-measurable, where B{R + ) is Borel a-algebra 

(ii) Cddldg cocycle property: <p(t,u) forms a (perfect) cddldg cocycle over 6 if it is 
cddldg in t and satisfies 

<p(0,cj) = idx, for all uj G Q, (1) 
<p(t + s,oj) = (p(t, 8 s u) o (p(s,u), (2) 
for all s, t G M+ and u G fl. 

A random dynamical system(RDS) induces a skew product flow of measurable maps 

9 t :0xX ^ flxX 

(uj,x) i — y (9 t cu,ip(t,u)xy 

The flow property Qt+s = Ot ° O s follows from ([2]). Denote by ^(Q x X) the probability 
measures on (Q x X, T®B). Moreover, Q t acts on fi G ^(llxX) by (9 t //)(C) = ^(O^C), 
for C G J 7 <g) B, t G R+. 

Definition 2.2. ^4 probability measure fi G ^(fi x X) called invariant for the skew 
product flow 6( if 

(i) the marginal of fi on Q is P, 

(ii) 0(/i = fi for all t G M + . 

If X is a Polish space with its Borel cr-algebra £>(X), every measure [i G &>(Q x X) with 
marginal P can be uniquely characterized through its factorization 

fi(du,dx) = /i^(dx)P(do;), 

where //^(dx) is a probability kernel, i.e. for any B G i3(X), fi.(B) is J-"-measurable; for 
P.a.s.u G fi, is a probability measure on (X, £>(X))(c.f.[3J p. 23]). Thus /i is invariant 
if and only if 

E^^c^/i.^J 7 ]^) = no t u), P.a.s, (3) 

for all t G M+. 

2.2. Sobolev spaces. Let Co(lR d ) be the space of continuous functions / on M. d satisfying 
lim f{x) = with norm H/H^rR*) = sup|/(x)|. Let C$(R d ) be the set of / G C (M. d ) 

\x\— ¥oo x 

such that / is 2 times different iable and the partial derivatives of / with order ^ 2 belong 
to C (M. d ). Let C™(R d ) stand for the space of all n times differentiable functions on M. 
with compact supports. Let >S(lR d ) be the Schwartz space of all rapidly decreasing real 
valued C°° functions on M d and S'(M. d ) the space of all tempered distributions on M d . Let 
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/ and / denote the Fourier transform and the Fourier inversion transform of / G S(R d ), 
respectively, i.e. 

for all u G R d . And Fourier transforms and the Fourier inversion transforms can be defined 
on S"(R d ) by the same means to the above one. We introduce the following Sobolev space 

HI A ' 2 (R d ) :={feS'(R d ): ||/|U, 2 < oo}, 

for any A G R, where 

\\f\\l^ m ■= [ (l + M 2 ) A |/(ii)| 2 du. 
In particular, H°- 2 (R d ) = L 2 (R d ). 
2.3. Rotation invariant a-stable process. 

Definition 2.3. A process L = (L t ) t>0 with L = a.s. is a d- dimensional Levy process 
if 

(i) L has independent increments; that independent of L v — L u if (u,v)C\ 
(s,*) = 0; 

(ii) L has stationary increments; that is, L t — L s has the same distribution as L v — L u 
if t — s = v — u > 0; 

(Hi) L t is stochastically continuous; 
(iv) L t is right continuous with left limit. 

Its characteristic function is given by 

E (exp{i(^, L t )}) = exp{*tf (z)}, z G R d . 

The function \1/ : R d — > C is called the characteristic exponent of the Levy process L. By 
Levy-Khintchine formula, there exist a nonnegative-definite d x d matrix Q, a measure v 
on R d satisfying 

z/({0}) = 0and / (|m| 2 A l)u(du) < oo, 

JR d \{0} 

and 7 G R d such that 

*(z) = i(zn) ~\(z,Qz) 

+ / {e^ -l-i(z,u)l {uHl )v(du). (4) 

JR d \{0} 

The measure v is called the Levy measure. 

Definition 2.4. For a G (0,2). A d-dimensional rotation invariant a-stable process L is 
a Levy process such that its characteristic exponent ^ is given by 

\fr(z) = -C\z\ a , zeR d . 

4 



Thus, for c?-dimensional rotation invariant a-stable process L, Levy measure v is given 

by 

and Q = in (J3J). Moreover, 



-C\z\* = I ( e ^>-l-i(^ w )l M<1 )_^d«. 



Define 



(C a f)(x) := / (/(X + M) - /(g) - (d x f(x),u)l\ u \^ 1 ) \i+a du 



on C^(M d ). And then for (Gl d 

(A^>)(x) = [ (e*<«*> - 1 - z-(£,«>lHa)r%^. 



By Courrege's second theorem (c.f.p] Theorem 3.5.5, p. 183]), for every / G C^°(lR a! ) 



(£»/)(*) = — Lj/ e^[e-^(£ a e^^)(x)]/»dz 



(2vr) d / 2 
1 



D i(z,x) 



(27f) d / 2 J Rd [J R d\ {0} " \u 

C 



f(z)dz 

\u\ 

= c.[-(-Ar /2 /]M. 

Moreover, the following result is well-known (c.f.jT]). 

Theorem 2.5. Let C a be as above for a G (0,2) and £2 = A, as defined on C'^'(]R d ) 
in L 2 (R d ). Then C a , < a ^ 2, has a unique closed extensions to self-adjoint negative 
operators on the domain H. a ' 2 (M. d ). 

3. From Markov measures to stationary measures 

Let (U, U, n) be a a-finite measurable space. Let {W(t)} t ^o be an m- dimensional stan- 
dard J-^-adapted Brownian motion, and {k t , t ^ 0} a stationary J-" r adapted Poisson point 
process with values in U and with characteristic measure n(cf. [7]). Let N k ((0, t], du) be 
the counting measure of k t , i.e., for A G U 

N k ((0,t],A) := #{0<s<:t:k s eA}, 

where # denotes the cardinality of a set. The compensator measure of N k is given by 

N k {(0, t],du) := JV fc ((0, t], du) - tn{du). 

Fix a Uq G U such that n(U — Uq) < 00, and consider the following SDE with jumps 



m S& d : 



X t {x) = x+ [ b{X s {x))ds + [ <j(X s (x))dW s 
Jo Jo 



+ ! + I f{X s 4x),u)N k (ds,du) 

JO JVo 

+ f + f g(X s 4x),u)N k (ds,du), t^O, (5) 

Jo JV-V 

where b : R d \-> R d } a : R d ^ R d x R m , f, g : R d x U m- R d satisfy the following 
assumptions: 

(Hfc) there exists a constant C& > such that for x, ?/ G R d 

\b(x) - b(y)\ < C 6 |a; - y\ ■ \og(\x - y^ 1 + e); 

(Ho-) there exists a constant C a > such that for x, y G K d 

- a(y)\ 2 ^ - y| 2 • logflar - y^ 1 + e); 

(Hj) for some g > (2d) V 4 and any p G [2, g], there exists a constant C p > such that 
for x,y eR d 

\f(x,u) - f(y,u)\ p n(du) ^ C p \x - y\ p ■ log(\x - y\~ l + e), 

'Do 

and 

/ \f(x,u)\ p n(du)^C p (l + \x\r. 

JVo 

(H ) for u G U — U , x i — y g(x, u) G C(R d ), where C(R d ) stands for the total of continuous 
functions from R d to R d . 

Here, the second integral of the right hand side in Eq.([5]) is taken in Ito's sense, and the 
definitions of the third and fourth integrals are referred to [7j. 

Under (lib), (Ho-), (H/) and (H 9 ), it is well known that there exists a unique strong 
solution to Eq.([5])(cf. [TO])- This solution will be denoted by X t (x). Set 

& >t := a{W s , N p ((0, s],B);s^t,Be U}, 
& <t := a{W s , N p ((0, s], B); s ^ t, B G U}, 

for t ^ 0. And then 
and 

@t C J^f+s, S #5 0. 

Moreover, X t (x) is J^-measurable and independent of 9^ l ^>o- 

Definition 3.1. A probability measure fi on (Q x R^J 7 ® ,B(R d )) is called a Markov 
measure if fi w satisfies 

E(//.|jF <QO ) =E(/i.|JT =0 ), P.a.s.. 
For the Markov process X t (x), the transition probability is defined by 
p t (x,B) := P(X t (x) eB), t>0,B G B(R d ). 

Definition 3.2. A measure ft on (R d , B(R d )) is a stationary measure for p t or Eq.^) if 
i p t {x, B)p(dx) = /2(B), Vt > 0, B G B(R d ). 

JR d 



Theorem 3.3. Set <p(t,u)x := X t (x), and then for an invariant Markov measure \i of 
the skew product flow Q t! the stationary measure p, for pt is given by 



/2 = E(//.|jr >0 ). 



Proof. By Definition 13.11 



Set 



E(/i.|#- <00 ) =E(//.|# 



* o) 



So), 



p, := E(//.|«^ ; 

and then it follows from independence of J^>o and J^=o 

fx = E[E(/i.|^ =0 )|^> ] 



E[E(/i.|^ =0 )] 
E(/i.)- 



Therefore v is not random. 
By it holds that for t > 



E[<£>(t, •)/^.|^ 1 #>o] = E\E((p(t, •)/x.|^ _1 J r )|^ 1 # > o] 

= E[ M .|^ 1 ^ > o] 

= E^.l^o)^) 

= /i. 



Besides, 



E^^-^.l^^o] = EfE^^-^.l^o)!^ 1 ^ 
= E^t,-)/^ 1 ^] 

Thus for G B(R d ) 

fc(B) = E[i f (t,-)fc(B)]=E[fc(if(t,-)- 1 B)] 



>oJ 



P(du) / l B (^(t,w)x)/i(dx) 

/i(da;) / l B ((p(t, u)x)P(dui) 
d Jq 

p t (x,B)fx(dx). 

d 

By Definition 13.21 E(/x, |J^n) is a stationary measure for □ 

So, this theorem, together with [101 Theorem 1.1] and [5], Lemma 5.1], yields 

Theorem 3.4. If (M. d ,B(M> d )) is replaced by its one point compactification (W d , B(W d )), 
stationary measures for Eq.$) on (R d ,B(R d )) exist. 
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4. From Fokker-Planck equations to stationary measures 
Consider the following equation 

X t (x)=x + [ b(X s (x))ds+ [ a(X s (x))dW s + L t , 
Jo Jo 

where L t is a rotation invariant a-stable process independent of W t . Based on the levy-Ito 
representation of L t , the above equation can be rewritten as follows: 

(i) for 1 ^ a < 2, 

X t {x) = x+ I b{X s {x))ds+ [ a(X s (x))dW s 
Jo Jo 

+ uN k (ds,du) + / / uN k (ds,du), (6) 

Jo J\u\^8 Jo J\u\>8 

where k(t) := L t — L t _ and < 5 < 1 satisfies that for some q > 4d 

q6 



(1 - 5)i+ l 

(ii) for < a < 1 



< 1; 



X t {x) = x+ I b(X s (x))ds+ / a{X s {x))dW s 



'0 ^0 

+ / / uN k (ds,du). (7) 



We study mainly Eq.([6]). Eq.Q can be dealt with similarly. 

Under (Hf,) and (H CT ), by jTOj, Theorem 1.3], for almost all u G ft, x \-> X t (x,u) is a 
homeomorphism mapping on Mr, where X t (x,u) is the solution of Eq.([6]). Define 

(p t h)(x):=E[h(X t (x))]= [ h{y)p t {x,dy), 

Jm. d 

for h G C (R d ). Thus p t h e C (R d ) by dominated convergence theorem. Let M r (M. d ) be 
the set of all finite regular signed measures on B{R d ). And then it is adjoint of Co(M d )(c.f. 
[HI Theorem 5.4.3, p. 148]. 

Lemma 4.1. The family of operators {pt,t ^ 0} defined above is a strongly continuous 
contraction semigroup on Co(K d ). 

Proof. For t, s > and h G C {R d ), by C-K equation 

ip t+s h)(x) = / h(y)p t +s(x,dy) = / h(y) / p s (z, dy)p t (x, dz) 



p t (x,dz) I h(y)p s (z,dy) = \ p t (x,dz)(p s h)(z) 
= (Pt(p s h))(x). 
So, pt+s = PtPs- Obviously p = I. 



Next, we prove strong continuity. For h G Co(ffi d ), h is uniformly continuous on M. d . 
And for We > 0, there exists an 77 > such that \h(x) —h(y)\ < e for x, y G M. d , \x — y\ < rj. 
For any A G M r {R d ) 



(p t h)(x)X(dx) — / h(x)X(dx 



h(y)p t (x,dy)X(dx) 



\h(y) - h{x)\p t (x,dy)\X\(dx) 



h(x)p t (x, dy)X(dx) 



\h(y) - h(x)\p t (x,dy)\X\(dx) 



x—y\<r] 



+ 



\h(y) - h(x)\p t (x,dy)\X\(dx) 



^ \X\(R d )e + 2\\h\\ / P{\X t (x)-x\^r]}\X\(dx), 

where |A| stands for the variation measure of the signed measure A. For J Rd P{\X t (x)—x\ ^ 
77} I A I (da:) , by stochastical continuity of X t and dominated convergence theorem, when t 
is small enough, 



P{\X t (x)-x\ ^ ?7}|A|(dx) < e. 



So, 



lim / (p t h)(x)X(dx) 
40 J R d 



h(x)X(dx) 



That is to say, p t h converges weakly to h. By [TBI Theorem, p. 233], p t h converges strongly 
to h. 

Finally, by [9j Definition 2.1, p. 4], {pt,t ^ 0} is a strongly continuous contraction 



semigroup on Co(R d ). 

Let £ be the infinitesimal generator of {pt,t ^ 0}. 
Lemma 4.2. For h G C 2 {R d ), 

1 / d 2 \ 
(£h)(x) = (d x h(x),b(x)) + - ^ — dx K x ) J °iA x ) + (£ a h)(x). 

Proof. By the Ito formula, one can obtain for h G C 2 (R d ) 



h(X t ) - h{x) 



□ 







(d y h(X s ),b(X s ))ds+ / (d y h(X s ),a(X s )dW s ) 







/ / (h(X s + u)-h(X s ))N p (ds,du) 

'0 J\u\45 

+ f [ (h(X s +u)- h(X s )) N p (ds, du) 

J0 J\u\>8 
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+ 



Hi 



o 2 



+ 




J\u\^S 



d y idyj h ^)^ dS 
(h(X s + u)- h(X s ) 



(d y h(X s ),u)) , d J ' a _ duds. 



\u 



d+a 



Taking expectation on two sides, we get 



(p t h)(x) - h(x) 



f E[(d y h(X s ),b(X s ))]ds + l f E 
o * Jo 

+ f e\ f (h(X s + u)-h(X s ) 



h(X s ) a tJ {X s ) 



■(dyh(X s ),U)) 



a 



d,a 



U 



d+a 



du 



ds 




JR d 



(d y h(y),b(y))p s (x,dy)ds 



+- 



+ 




d 2 



Kv) °ij(.y)Ps(x,dy)ds 



(%<% ' 

(h(y + u) - h(y) 
C 

-(d y h(y),u)) t A d d * a du p s (x,dy)ds 




JR d 



\U 



Thus 



1 



lim- ((p t h)(x) -h{x)) = (d x h(x), b(x)) + - 
40 t 2 \ oxiOXj 



(J 2 



h(x) CTij(x) 



+ 



M d \{0} 



(h(x + u) — h(x) — (d x h(x) , u)) 



C A 



\u 



d+a 



du. 



By [IH Lemma 31.7, p.209] 

(Ch)(x) = (d x h{.r).bi.r)}-- 

+ 



d 2 



2 \ dxidxj 



h(x) aij(x) 



(h(x + u) - h(x) - (d x h(x),u)) '° Q du 



(d x h(x),b(x)} + - 



1 / d 2 



2 \dxidxj 



h(x) <Jij(x) + (C a h)(x). 



To get the main result, we make the following assumptions. 
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(Hp) For all t > and x G M d , the transition probability p t (x,dy) admits a density 
p t (x,y) and functions 

(t,v) ^ ^Pt(x,y), 
d 

(t,y) ^ Q—(bj(y)p t (x,y)),j = 1,2, 

(*>!/) H> g^g^-(<T ij (y)pt(x 1 y)),i 1 j = 1,2, . . . ,d, 

i-). (C a p t (x, -))(y), 

exist and are continuous on M + x M d . 
Under (H p ), the distribution of X t (x) has a density px(t,y) = p t (x,y) for £ > 0. 

Lemma 4.3. px(t,y) satisfies the following Fokker-Planck equation 

jtu(t,y) = (jr*u(t,-))(y), 

lim«(t, y)dy = 4(dy), I 8 ) 
where C* is adjoint of C, i.e. 

for<pe C c °°(M d ). 

Proof. By similar deduction to that of Lemma [4.21 we get for ^ e C£°(K d ) and any r > 

Hv)Px{t + r,y)dy - / ip(y)p x (t,y)dy 
d Jm d 

t+r 



(d y ip(y), b(y))p x {s, y)dyds 

rt+r p 

+ / / {£a^){y)px{s,y)dyds. 
Jt J«. d 

Divided by r and Letting r 1 0, it follows from dominated convergence theorem and 
integration by parts that 

d f 

i>(y)-^-px(t,y)dy = / (d y ip(y),b(y))p x (t,y)dy 

at JR d 

+1 2LX^k m ) a,iiy)Mt ' y)Ay 

+ / (C a i))(y)px(t,y)dy 

JR d 

= ~ / ip{y)w-(bj{y)p x {t,y))dy 
jR d oy,j 
11 



+ / ^(y)(C aPx (t,-))(y)dy. 



Thus 

d d Id 2 

foPx(t,y) = --7r-(bj(y)px(t,y)) + 2dy~~dy~^ 3 ^ Px ^ ,y ^ + 0)(v) 

= (£*p x (t,-))(y)- 



Note that Schertzer et. al. [12] also derived the Fokker-Planck equation in this context. 
The main result in this section is the following theorem. 

Theorem 4.4. If p(y) G H 2 ' 2 (IR d ) satisfies the following equation 



Then p(dy) := p(y)dy is a stationary measure for pt. 

Proof. By Lemma [4.11 and [4.21 {pt,t ^ 0} is a strongly continuous semigroup on Co(M. d ) 
with the infinitesimal generator C. Let pi be adjoint of pt- By Theorem I2.5[ H 2 ' 2 (R d ) 
is the closure of V(£*) in M r (M d )- [HI Theorem 10.4, p. 41] admits us to get that the 
restriction pf of p* to H 2 ' 2 (M d ) is a strongly continuous semigroup on M r (M. d ). Moreover, 
the infinitesimal generator C + of pf is the part of C* in HI 2 ' 2 (lR d ), i.e. T>{£, + ) = {h G 
£>(£*) n M 2 ' 2 (R d ),£*h G U 2 ' 2 (R d )} and C+h = C*h for h G £>(£+). Thus, Eq.® has a 
unique solution in H 2 ' 2 (M d ) by Theorem 1.3, p.102]. 
Next, since p(y) satisfies 



□ 



C*p = 



and 





C*p = 



and 





p(y) solves Eq.® and is a density function. So, 

Px(t,y) = p(y), t>0. 

For t > and B G 
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px(s + t,y)dy 

By Definition 13.21 p(dy) is a stationary measure for p t . □ 



B 



B 



Remark 4.5. If b(x) = —x and c(x) = 0, the above theorem is [U Proposition 3.2(H)]. 

Moreover, 

p(u) = exp{--C>| Q }, u G R d , 
a 

where C is the same constant as one in Definition \2.4 



5. From probability density functions to stationary measures 

In the section we study Eq.(j5]) with g = 0. We further make the following assumptions. 
(H^ CT j) b and a are (4c? + 6)-times differentiable with bounded derivatives of all order 
between 1 and 4g? + 6. Besides, f{-,u) is (4d + 6)-times differentiable, and 

/((),-) G p| L«(U ,n) 

2^<?<oo 

swp\d r J(x,-)\ G p| L q (U Q ,n), l^r^4rf + 6, 

X 2sg<j<°o 

where the space (Uo,n) is equipped with a norm and d r x f(x, •) stands for r order partial 
derivative of f(x, •) with respect to x. 

(H^ CT j) There exist three constants e > 0, 5 ^ and C > such that for all x, y G R d 

(y, cr(x)o- T (x)y) ^ |y 

and 



1 + |x| 

det{/ + ra x /(x,M)}| >C 



for all r G [0, 1]. 

Under (H^j) and (H| by [U Theorem 2-29, p. 15], Eq.([5]) has a unique solution 
denoted by X t . Moreover, the transition probability p t (x,dy) has a density p t (x,y) and 
(t,x,y) i — ^ p t (x,y) is continuous. Thus, the distribution of JY" t has a density px(t,y) = 
Pt(x,y). 

Theorem 5.1. Suppose that lim p x (t,y) = p(y), where p(y) satisfies 

t—>oo 

p(y)^0,\/yeR d and [ p{y)dy = 1. 

T/ien p{dy) := p(y)dy is a stationary measure for p t . 
Proof. For t > and £> G i3(R d ), 

p t (x, B)p(dx) = / pt(x, i?)p(x)dx = / p t (x, B) lim px(s, x)dx 

jRd J R d s^oo 

= lim / p t (x, B)p x (s,x)dx = lim / / p t (x, y)px(s, x)dydx 
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= lim / p x (s + t,y)dy 



= lim / / p t (x,y)px(s,x)dxdy 




= KB). 



By Definition 13.21 p(dy) is a stationary measure for p t . 



□ 



Remark 5.2. By the above theorem, we see that if a limiting distribution exists, it must 
be a stationary measure. This theorem also has a corresponding version in the theory of 
Markov chains (c.f. [8, p. 237],). 
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